Question 2: Focus on the evolution of SDG scores over time

data_question2 <- read.csv(here("scripts", "data", "data_question24.csv"))
data_question2 <- data_question2 %>% select(-X)

General time evolution of SDG socres

First, we look at the evolution of SDG achievement overall score over time for the whole world, by continent and by region.

data1 <- data_question2 %>% group_by(year) %>%
  mutate(mean_overall_score_by_year=mean(overallscore))

ggplot(data1) +
  geom_line(mapping=aes(x=year, y=mean_overall_score_by_year), color="blue", lwd=0.8) +
  geom_point(mapping=aes(x=year, y=mean_overall_score_by_year), color="blue", lwd=2.5) +
  scale_y_continuous(limits = c(0, 100))  +
  labs(title = "Slow evolution of overall SDG scores across the world",
       y = "Mean Overall SDG Score",
       x = "Year"
       )

The general evolution of SDG scores around the world is increasing over the years, but very slowly.

data2 <- data_question2 %>% group_by(year, continent) %>%
  mutate(mean_overall_score_by_year=mean(overallscore))

ggplot(data2) +
  geom_line(mapping=aes(x=year, y=mean_overall_score_by_year, color=continent), lwd=0.8) +
  geom_point(mapping=aes(x=year, y=mean_overall_score_by_year, color=continent), lwd=1.5) +
  scale_y_continuous(limits = c(0, 100)) +
  labs(title = "Evolution of the mean overall SDG achievement score by continent",
       y = "Mean Overall SDG Score",
       x = "Year"
       )

Looking at the continents, we see that Europe is above the others, while Africa is below, but in general, all have increasing overall scores.

data3 <- data_question2 %>% group_by(year, region) %>%
  mutate(mean_overall_score_by_year=mean(overallscore))

ggplot(data3) +
  geom_line(mapping=aes(x=year, y=mean_overall_score_by_year, color=region), lwd=0.8) +
  geom_point(mapping=aes(x=year, y=mean_overall_score_by_year, color=region), lwd=1.5) +
  scale_y_continuous(limits = c(0, 100)) +
  labs(title = "Evolution of the mean overall SDG achievement score by region",
       y = "Mean Overall SDG Score",
       x = "Year"
       )

This view that groups the countries by region gives us precision about the previous information. Indeed, it is Western Europe that is particularly above and Sub-Saharan Africa that is clearly below.

Second, we look at the evolution of SDG achievement scores(16) over time for the whole world, by continent and by region.

data4 <- data_question2 %>%
  group_by(year) %>%
  summarise(across(starts_with("goal"), mean, na.rm=TRUE)) %>%
  pivot_longer(cols = starts_with("goal"), names_to = "goal", values_to = "mean_value")

color_palette <- c("red", "blue", "green", "orange", "purple", "pink", "lightblue", "gray", "cyan", "magenta", "yellow", "darkgreen", "darkblue", "darkred", "darkgrey", "darkcyan")

ggplot(data = data4) +
  geom_line(mapping = aes(x = year, y = mean_value, color = goal), size = 0.7) +
  geom_point(mapping = aes(x = year, y = mean_value, color = goal), size = 1) +
  scale_color_manual(values = color_palette) +
  scale_y_continuous(limits = c(0, 100)) +
  labs(title = "Evolution of the mean SDG achievement scores across the world",
       y = "Mean SDG Scores",
       x = "Year"
       ) +
  guides(
    color = guide_legend(
      ncol = 2,        # Number of columns
      title.position = "top",  # Position of the legend title
      title.hjust = 0.5        # Horizontal justification of the legend title
    )
  )

Here, by looking at the SDGs individually, we notice that all SDGs except from goal 9 (industry innovation and infrastructure) are close to one another in terms of level and growth. Goal 9 starts far below the others in 2000 and growths faster until almost exceeding 50%.

color_palette <- c("red", "blue", "green", "orange", "purple", "pink", "lightblue", "gray", "cyan", "magenta", "yellow", "darkgreen", "darkblue", "darkred", "darkgrey", "darkcyan")

ggplot(data = data4) +
  geom_line(mapping = aes(x = year, y = mean_value, color=goal), size = 0.8) +
  geom_point(mapping = aes(x = year, y = mean_value, color=goal), size = 1.5) +
  scale_color_manual(values = color_palette) +
  scale_y_continuous(limits = c(0, 100)) +
  labs(title = "Evolution of the mean SDG achievement scores across the world",
       y = "Mean SDG Scores",
       x = "Year"
       ) +
  facet_wrap(~ goal, nrow = 4)

In contrast to the aspect discussed in the precedent graph, this graph shows us the same information in a different way and it pops out that some goals did not increase their scores much in the last two decades, for example goal 13 (climate action) and goal 12 (responsible consumption and production).

data5 <- data_question2 %>%
  group_by(year, continent) %>%
  summarise(across(starts_with("goal"), mean, na.rm=TRUE)) %>%
  pivot_longer(cols = starts_with("goal"), names_to = "goal", values_to = "mean_value")

color_palette <- c("red", "blue", "green", "orange", "purple", "pink", "brown", "gray", "cyan", "magenta", "yellow", "darkgreen", "darkblue", "black", "darkorange", "darkcyan")

ggplot(data = data5) +
  geom_line(mapping = aes(x = year, y = mean_value, color=continent), size = 0.7) +
  scale_color_manual(values = color_palette) +
  scale_y_continuous(limits = c(0, 100)) +
  labs(title = "Evolution of the mean SDG achievement scores by continent",
       y = "Mean SDG Scores",
       x = "Years from 2000 to 2022"
       ) +
  facet_wrap(~ goal, nrow = 4)+
  scale_x_continuous(breaks = NULL)+
  theme_light()

Now, comparing the SDG scores by continent, we observe that most of the time, Europe is at the top of the graph and Africa at the bottom, exept for goals 12 and 13 that are linked to ecology. Some other information stand out:

  • Americas are far behind the other parts of the world regarding goal 10: reduced inequalities.

  • Africa is far behind the other continents (even if becoming better) for goals 1, 3, 4 and 7.

  • Goal 9 (industry, innovation and infrastructure) show exponential growth for almost all continents.

Third we create an interactive map of the world to be able to navigate from year 2000 to 2022, seeing the level of achievement of the SDGs (overall score) for each country. To open it in your browser, use this R file: interactive_map_1. Here is only a non-interactive world map of the overall SDGs achievement scores, not taking into account the evolution over the years.

library(rnaturalearth)
library(tidyverse)
library(sf)
# Load world map data
world <- ne_countries(scale = "medium", returnclass = "sf")

# Merge data with the world map data
data0 <- merge(world, data_question2, by.x = "iso_a3", by.y = "code", all.x = TRUE)

data0 %>%
  sf::st_transform(crs = "+proj=robin") %>%
  ggplot() +
  geom_sf(color = "lightgrey") +
  geom_sf(aes(fill = overallscore), color = NA) +
  scale_fill_gradientn(
    colors = c("darkred", "orange", "yellow", "darkgreen"),
    values = scales::rescale(c(0, 0.25, 0.5, 1)),
    name = "Overall Score",
    na.value = NA
  ) +
  labs(title = "Mean overall SDG achievement score by country")+
  coord_sf(datum = NA) +
  theme_minimal()

SDG adoption in 2015

binary2015 <- data_question2 %>% 
  group_by(code) %>%
  mutate(across(5:21, ~ . - dplyr::lag(.), .names = "diff_{.col}")) %>%
  ungroup()

# Create a new column (binary variable) with value 1 if the year is after 2015 and zero otherwise. 
binary2015 <- binary2015 %>% 
  mutate(after2015 = ifelse(year > 2015, 1, 0)) %>%
  filter(as.numeric(year)>=2009)

Look at the distribution of the difference in SDG scores from one year to the next.

# histogram of difference in scores between years
unique_years <- unique(binary2015$year)
library(viridis)
plotly::plot_ly() %>%
  plotly::add_trace(
    type = "histogram", 
    data = binary2015, 
    x = ~diff_overallscore,
    color = ~diff_overallscore
  ) %>%
  plotly::layout(
    title = "Distribution of SDG evolution",
    xaxis = list(title = "Year difference SDG score", range = c(-3,3)),
    yaxis = list(title = "Frequency", range = c(0,40)),
    sliders = list(
      list(
        active = 0,
        currentvalue = list(prefix = "Year: "),
        steps = lapply(seq_along(unique_years), function(i) {
          year <- unique_years[i]
          list(
            label = as.character(year),
            method = "restyle",
            args = list(
              list(x = list(binary2015$diff_overallscore[binary2015$year == year]))
            )
          )
        })
      )
    )
  )

2 fitted regression lines (before and after 2015) to see if there is a jump after the adoption and if the the SDG scores increase faster.

# Graphs to show the jump (or not) in 2015

# Filter data
data_after_2015 <- filter(binary2015, as.numeric(year) >= 2015)
data_before_2016 <- filter(binary2015, as.numeric(year) <= 2015)

plotly::plot_ly() %>%
  plotly::add_trace(data = data_after_2015, x = ~year, y = ~fitted(lm(overallscore ~ year, data = data_after_2015)), type = 'scatter', mode = 'lines', line = list(color = 'blue'), name = "After 2015") %>%
  plotly::add_trace(data = data_before_2016, x = ~year, y = ~fitted(lm(overallscore ~ year, data = data_before_2016)), type = 'scatter', mode = 'lines', line = list(color = 'red'), name = "Before 2015") %>%
  plotly::layout(title = "Different patterns across SDGs before and after 2015",
         xaxis = list(title = "Year"),
         yaxis = list(title = "SDG achievement score", range = c(30, 85)),
         shapes = list(
           list(
             type = 'line',
             x0 = 2015,
             x1 = 2015,
             y0 = 0,
             y1 = 1,
             yref = 'paper',
             line = list(color = 'grey', width = 2, dash = 'dot')
           )
         ),
         updatemenus = list(
           list(
             buttons = list(
               list(
                 args = list("y", list(
                   ~fitted(lm(overallscore ~ year, data = data_after_2015)),
                   ~fitted(lm(overallscore ~ year, data = data_before_2016))
                 )),
                 label = "Overall score",
                 method = "restyle"
               ),
               list(
                 args = list("y", list(
                   ~fitted(lm(goal1 ~ year, data = data_after_2015)),
                   ~fitted(lm(goal1 ~ year, data = data_before_2016))
                 )),
                 label = "Goal 1",
                 method = "restyle"
               ),
               list(
                 args = list("y", list(
                   ~fitted(lm(goal2 ~ year, data = data_after_2015)),
                   ~fitted(lm(goal2 ~ year, data = data_before_2016))
                 )),
                 label = "Goal 2",
                 method = "restyle"
               ),
               list(
                 args = list("y", list(
                   ~fitted(lm(goal3 ~ year, data = data_after_2015)),
                   ~fitted(lm(goal3 ~ year, data = data_before_2016))
                 )),
                 label = "Goal 3",
                 method = "restyle"
               ),
               list(
                 args = list("y", list(
                   ~fitted(lm(goal4 ~ year, data = data_after_2015)),
                   ~fitted(lm(goal4 ~ year, data = data_before_2016))
                 )),
                 label = "Goal 4",
                 method = "restyle"
               ),
               list(
                 args = list("y", list(
                   ~fitted(lm(goal5 ~ year, data = data_after_2015)),
                   ~fitted(lm(goal5 ~ year, data = data_before_2016))
                 )),
                 label = "Goal 5",
                 method = "restyle"
               ), 
               list(
                 args = list("y", list(
                   ~fitted(lm(goal6 ~ year, data = data_after_2015)),
                   ~fitted(lm(goal6 ~ year, data = data_before_2016))
                 )),
                 label = "Goal 6",
                 method = "restyle"
               ),
               list(
                 args = list("y", list(
                   ~fitted(lm(goal7 ~ year, data = data_after_2015)),
                   ~fitted(lm(goal7 ~ year, data = data_before_2016))
                 )),
                 label = "Goal 7",
                 method = "restyle"
               ),
               list(
                 args = list("y", list(
                   ~fitted(lm(goal8 ~ year, data = data_after_2015)),
                   ~fitted(lm(goal8 ~ year, data = data_before_2016))
                 )),
                 label = "Goal 8",
                 method = "restyle"
               ),
               list(
                 args = list("y", list(
                   ~fitted(lm(goal9 ~ year, data = data_after_2015)),
                   ~fitted(lm(goal9 ~ year, data = data_before_2016))
                 )),
                 label = "Goal 9",
                 method = "restyle"
               ), 
               list(
                 args = list("y", list(
                   ~fitted(lm(goal10 ~ year, data = data_after_2015)),
                   ~fitted(lm(goal10 ~ year, data = data_before_2016))
                 )),
                 label = "Goal 10",
                 method = "restyle"
               ),
               list(
                 args = list("y", list(
                   ~fitted(lm(goal11 ~ year, data = data_after_2015)),
                   ~fitted(lm(goal11 ~ year, data = data_before_2016))
                 )),
                 label = "Goal 11",
                 method = "restyle"
               ),
               list(
                 args = list("y", list(
                   ~fitted(lm(goal12 ~ year, data = data_after_2015)),
                   ~fitted(lm(goal12 ~ year, data = data_before_2016))
                 )),
                 label = "Goal 12",
                 method = "restyle"
               ),
               list(
                 args = list("y", list(
                   ~fitted(lm(goal13 ~ year, data = data_after_2015)),
                   ~fitted(lm(goal13 ~ year, data = data_before_2016))
                 )),
                 label = "Goal 13",
                 method = "restyle"
               ), 
               list(
                 args = list("y", list(
                   ~fitted(lm(goal15 ~ year, data = data_after_2015)),
                   ~fitted(lm(goal15 ~ year, data = data_before_2016))
                 )),
                 label = "Goal 15",
                 method = "restyle"
               ),
               list(
                 args = list("y", list(
                   ~fitted(lm(goal16 ~ year, data = data_after_2015)),
                   ~fitted(lm(goal16 ~ year, data = data_before_2016))
                 )),
                 label = "Goal 16",
                 method = "restyle"
               ),
               list(
                 args = list("y", list(
                   ~fitted(lm(goal17 ~ year, data = data_after_2015)),
                   ~fitted(lm(goal17 ~ year, data = data_before_2016))
                 )),
                 label = "Goal 17",
                 method = "restyle"
               )
             )
           )
         )
  )

Simple OLS regression on the difference between years of SDG scores:

# Simple linear regression of the overall score on the difference in SDG scores variables "after2015"
reg2.1 <- lm(diff_overallscore ~ after2015, data=binary2015)
reg2.1.1 <- lm(diff_goal1 ~ after2015, data=binary2015)
reg2.1.2 <- lm(diff_goal2 ~ after2015, data=binary2015)
reg2.1.3 <- lm(diff_goal3 ~ after2015, data=binary2015)
reg2.1.4 <- lm(diff_goal4 ~ after2015, data=binary2015)
reg2.1.5 <- lm(diff_goal5 ~ after2015, data=binary2015)
reg2.1.6 <- lm(diff_goal6 ~ after2015, data=binary2015)
reg2.1.7 <- lm(diff_goal7 ~ after2015, data=binary2015)
reg2.1.8 <- lm(diff_goal8 ~ after2015, data=binary2015)
reg2.1.9 <- lm(diff_goal9 ~ after2015, data=binary2015)
reg2.1.10 <- lm(diff_goal10 ~ after2015, data=binary2015)
reg2.1.11 <- lm(diff_goal11 ~ after2015, data=binary2015)
reg2.1.12 <- lm(diff_goal12 ~ after2015, data=binary2015)
reg2.1.13 <- lm(diff_goal13 ~ after2015, data=binary2015)
reg2.1.15 <- lm(diff_goal15 ~ after2015, data=binary2015)
reg2.1.16 <- lm(diff_goal16 ~ after2015, data=binary2015)
reg2.1.17 <- lm(diff_goal17 ~ after2015, data=binary2015)

stargazer::stargazer(reg2.1, reg2.1.1, reg2.1.2, reg2.1.3, reg2.1.4, reg2.1.5,
          title="Impact of the 2015 adoption of SDG by the UN (Goals 1-5)",
          type='html',
          column.labels=c("Overall score", "Goal 1", "Goal2", "Goal 3", "Goal 4","Goal 5"),
          digits=3,
          omit = "coef",
          dep.var.labels.include = FALSE,
          model.numbers = FALSE

)
Impact of the 2015 adoption of SDG by the UN (Goals 1-5)
Dependent variable:
Overall score Goal 1 Goal2 Goal 3 Goal 4 Goal 5
after2015 -0.073*** -0.194** -0.116 -0.423*** -0.275** -0.315***
(0.024) (0.091) (0.072) (0.050) (0.125) (0.077)
Constant 0.398*** 0.535*** 0.259*** 0.764*** 0.582*** 0.863***
(0.017) (0.065) (0.051) (0.035) (0.088) (0.055)
Observations 2,170 2,002 2,170 2,170 2,170 2,170
R2 0.004 0.002 0.001 0.033 0.002 0.008
Adjusted R2 0.004 0.002 0.001 0.032 0.002 0.007
Residual Std. Error 0.570 (df = 2168) 2.042 (df = 2000) 1.685 (df = 2168) 1.154 (df = 2168) 2.900 (df = 2168) 1.803 (df = 2168)
F Statistic 8.957*** (df = 1; 2168) 4.496** (df = 1; 2000) 2.570 (df = 1; 2168) 72.875*** (df = 1; 2168) 4.892** (df = 1; 2168) 16.587*** (df = 1; 2168)
Note: p<0.1; p<0.05; p<0.01
stargazer::stargazer(reg2.1.6, reg2.1.7, reg2.1.8, reg2.1.9, reg2.1.10, reg2.1.11,
          title = "Impact of the 2015 adoption of SDG by the UN (Goals 6-11)",
          type = 'html',
          column.labels = c("Goal 6", "Goal 7", "Goal 8", "Goal 9", "Goal 10", "Goal 11"),
          digits = 3,
          omit = "coef",
          dep.var.labels.include = FALSE,
          model.numbers = FALSE)
Impact of the 2015 adoption of SDG by the UN (Goals 6-11)
Dependent variable:
Goal 6 Goal 7 Goal 8 Goal 9 Goal 10 Goal 11
after2015 -0.127*** -0.259*** 0.203*** -0.010 -0.130 -0.044
(0.018) (0.068) (0.048) (0.096) (0.187) (0.116)
Constant 0.276*** 0.510*** 0.072** 1.496*** 0.491*** 0.244***
(0.013) (0.048) (0.034) (0.068) (0.132) (0.082)
Observations 2,170 2,170 2,170 2,170 2,002 2,170
R2 0.021 0.007 0.008 0.00001 0.0002 0.0001
Adjusted R2 0.021 0.006 0.008 -0.0005 -0.0003 -0.0004
Residual Std. Error 0.428 (df = 2168) 1.588 (df = 2168) 1.112 (df = 2168) 2.246 (df = 2168) 4.179 (df = 2000) 2.700 (df = 2168)
F Statistic 47.616*** (df = 1; 2168) 14.432*** (df = 1; 2168) 18.066*** (df = 1; 2168) 0.011 (df = 1; 2168) 0.483 (df = 1; 2000) 0.146 (df = 1; 2168)
Note: p<0.1; p<0.05; p<0.01
stargazer::stargazer(reg2.1.12, reg2.1.13, reg2.1.15, reg2.1.16, reg2.1.17,
          title = "Impact of the 2015 adoption of SDG by the UN (Goals 12-17)",
          type = 'html',
          column.labels = c("Goal 12", "Goal 13", "Goal 15", "Goal 16", "Goal 17"),
          digits = 3,
          omit = "coef",
          dep.var.labels.include = FALSE,
          model.numbers = FALSE)
Impact of the 2015 adoption of SDG by the UN (Goals 12-17)
Dependent variable:
Goal 12 Goal 13 Goal 15 Goal 16 Goal 17
after2015 -0.002 -0.014 -0.225*** -0.260*** 0.709***
(0.034) (0.057) (0.063) (0.068) (0.088)
Constant 0.053** 0.094** 0.207*** 0.063 0.107*
(0.024) (0.041) (0.044) (0.048) (0.062)
Observations 2,170 2,170 2,170 2,170 2,170
R2 0.00000 0.00003 0.006 0.007 0.029
Adjusted R2 -0.0005 -0.0004 0.005 0.006 0.029
Residual Std. Error (df = 2168) 0.781 1.339 1.463 1.585 2.040
F Statistic (df = 1; 2168) 0.006 0.057 12.769*** 14.656*** 65.560***
Note: p<0.1; p<0.05; p<0.01

DiD using panel data:

# Create a panel data object
panel_data <- plm::pdata.frame(binary2015, index = c("country", "year"))

# Run the difference-in-differences model to take into account the general evolution over the years
reg2.2 <- plm::plm(diff_overallscore ~ after2015 + year + after2015:year, 
                 data = panel_data,
                 model = "within")
reg2.2.1 <- plm::plm(diff_goal1 ~ after2015 + year + after2015:year, 
              data = panel_data,
              model = "within")
reg2.2.2 <- plm::plm(diff_goal2 ~ after2015 + year + after2015:year, 
                data = panel_data,
                model = "within")
reg2.2.3 <- plm::plm(diff_goal3 ~ after2015 + year + after2015:year, 
                data = panel_data,
                model = "within")
reg2.2.4 <- plm::plm(diff_goal4 ~ after2015 + year + after2015:year, 
                data = panel_data,
                model = "within")
reg2.2.5 <- plm::plm(diff_goal5 ~ after2015 + year + after2015:year, 
                data = panel_data,
                model = "within")
reg2.2.6 <- plm::plm(diff_goal6 ~ after2015 + year + after2015:year, 
                data = panel_data,
                model = "within")
reg2.2.7 <- plm::plm(diff_goal7 ~ after2015 + year + after2015:year, 
                data = panel_data,
                model = "within")
reg2.2.8 <- plm::plm(diff_goal8 ~ after2015 + year + after2015:year, 
                data = panel_data,
                model = "within")
reg2.2.9 <- plm::plm(diff_goal9 ~ after2015 + year + after2015:year, 
                data = panel_data,
                model = "within")
reg2.2.10 <- plm::plm(diff_goal10 ~ after2015 + year + after2015:year, 
                data = panel_data,
                model = "within")
reg2.2.11 <- plm::plm(diff_goal11 ~ after2015 + year + after2015:year, 
                data = panel_data,
                model = "within")
reg2.2.12 <- plm::plm(diff_goal12 ~ after2015 + year + after2015:year, 
                data = panel_data,
                model = "within")
reg2.2.13 <- plm::plm(diff_goal13 ~ after2015 + year + after2015:year, 
                data = panel_data,
                model = "within")
reg2.2.15 <- plm::plm(diff_goal15 ~ after2015 + year + after2015:year, 
                data = panel_data,
                model = "within")
reg2.2.16 <- plm::plm(diff_goal16 ~ after2015 + year + after2015:year, 
                data = panel_data,
                model = "within")
reg2.2.17 <- plm::plm(diff_goal17 ~ after2015 + year + after2015:year, 
                data = panel_data,
                model = "within")

stargazer::stargazer(reg2.2, reg2.2.1, reg2.2.2, reg2.2.3, reg2.2.4, reg2.2.5,
          title="Impact of the 2015 adoption of SDG by the UN (Goals 1-5)",
          type='html',
          column.labels=c("Overall score", "Goal 1", "Goal2", "Goal 3", "Goal 4","Goal 5"),
          digits=3,
          omit = "coef",
          dep.var.labels.include = FALSE,
          model.numbers = FALSE
          
)
## 
## % Error: Unrecognized object type.
## % Error: Unrecognized object type.
## % Error: Unrecognized object type.
## % Error: Unrecognized object type.
## % Error: Unrecognized object type.
## % Error: Unrecognized object type.
stargazer::stargazer(reg2.2.6, reg2.2.7, reg2.2.8, reg2.2.9, reg2.2.10, reg2.2.11,
          title = "Impact of the 2015 adoption of SDG by the UN (Goals 6-11)",
          type = 'html',
          column.labels = c("Goal 6", "Goal 7", "Goal 8", "Goal 9", "Goal 10", "Goal 11"),
          digits = 3,
          omit = "coef",
          dep.var.labels.include = FALSE,
          model.numbers = FALSE)
## 
## % Error: Unrecognized object type.
## % Error: Unrecognized object type.
## % Error: Unrecognized object type.
## % Error: Unrecognized object type.
## % Error: Unrecognized object type.
## % Error: Unrecognized object type.
stargazer::stargazer(reg2.2.12, reg2.2.13, reg2.2.15, reg2.2.16, reg2.2.17,
          title = "Impact of the 2015 adoption of SDG by the UN (Goals 12-17)",
          type = 'html',
          column.labels = c("Goal 12", "Goal 13", "Goal 15", "Goal 16", "Goal 17"),
          digits = 3,
          omit = "coef",
          dep.var.labels.include = FALSE,
          model.numbers = FALSE)
## 
## % Error: Unrecognized object type.
## % Error: Unrecognized object type.
## % Error: Unrecognized object type.
## % Error: Unrecognized object type.
## % Error: Unrecognized object type.